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We show that there are effective three- and higher-body interactions generated by the two-body collisions of 
atoms confined in the lowest vibrational states of a 3D optical lattice. The collapse and revival dynamics of 
approximate coherent states loaded into a lattice are a particularly sensitive probe of these higher-body interac- 
tions; the visibility of interference fringes depend on both two-, three-, and higher-body energy scales, and these 
produce an initial dephasing that can help explain the surprisingly rapid decay of revivals seen in experiments. 
If inhomogeneities in the lattice system are sufficiently reduced, longer timescale partial and nearly full revivals 
will be visible. Using Feshbach resonances or control of the lattice potential it is possible to tune the effective 
higher-body interactions and simulate effective field theories in optical lattices. 
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I. INTRODUCTION 

The collapse and revival of matter-wave coherence is an expected consequence of two-body atom-atom interactions in trapped 
Bose-Einstein condensates (BECs) ||T1|2]|3]|4 |. Collapse and revival of few-atom coherent states in optical lattices has been seen 
in a number of experiments, first in single-well lattices |5 1 and subsequently in double- well lattices [|6||7]. In these experiments, 
a BEC is quickly loaded into a fairly deep 3D lattice such that the quantum state approximately factors into a product of coherent 
states localized to each lattice site ||8l|9]|T0l. Each coherent state, which is a superposition of different atom-number states, 
initially has a well-defined phase. If the lattice potential is quickly turned off before atom-atom interactions have a significant 
influence, the coherent states released from confinement at each site expand and overlap resulting in interference fringes in the 
imaged atom-density. However, if the atoms are held in the lattice for a longer duration before release, interactions will play a 
significant role by causing the phases of the different atom-number states in the superposition at each site to evolve at different 
rates. This will result in a dephasing of the coherent state, and a subsequent collapse of the interference fringe visibility after the 
atoms are released. For atoms in a homogenous lattice with two-body interactions and negligible tunneling, the coherent states 
at each lattice site are predicted to revive when the atom-number component states simultaneously re-phase after multiples of 
the time t2 = 2Trh/U2, where U2 is the two-body interaction energy H] El |3] IH |5| . 

In addition to the expected two-body physics described above, we show that the data in |5 , 6, 7| should also contain strong 
signatures of coherent three- and higher-body interactions. In contrast to the coherent dynamics described in this paper, re- 
cent experiments have studied inelastic three-body processes, including recent observations of Efimov physics ifTTl fT2l [131 . 
by tracking atom loss from recombination flTl. There has been a growing interest in three- and four-body physics (e.g., 
ifTSi [T6l [TTl [TSl [T9I I20I ). and the role of intrinsic three-body interactions on equilibrium quantum phases in optical lattices 
has been studied in 11211 122ll23l . The influence of higher bands on the Mott-insulator phase transition has been analyzed in ll24l . 
and three-body interactions of fermions and polar molecules in lattices have also been explored f25l. 

In this paper, we use the ideas of effective field theory to show that virtual transitions to higher vibrational states generate 
effective, coherent three-body interactions between atoms in the lowest vibrational states of a deep 3D lattice where tunneling 
can be neglected. More generally, virtual excitations also generate effective four- and higher-body interactions giving the non- 
equilibrium dynamics multiple energy scales. We also show that loading coherent states into an optical lattice creates a sensitive 
interferometer for probing higher-body interactions. In a sufficiently uniform lattice, multiple frequencies manifested as beatings 
in the visibility of the collapse and revival oscillations give a direct method for measuring the energy and frequency scales 
for elastic higher-body interactions. Remarkably, multiple-frequency collapse and revival patterns have been seen in recent 
experiments \2S\. 

Three-body interactions can also explain the surprisingly rapid damping of revivals seen in |j5]|6]|2l, where the overall visibility 
of the interference fringes decays after roughly 5 revivals (~ 3 ms for the system parameters in 016] |2l). This short timescale 
cannot be explained in terms of tunneling or atom loss. For example, for the system parameters in 15] IS I?), the tunneling- 
induced decoherence timescale has been found to be a factor of 10-100 times too long 1.27,1 . and the atom loss from three-body 
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recombination f]A\ appears to be negligible ll26ll . The latter observation is consistent with the expected three-body recombination 
timescales for *^Rb in a lattice [28 1. 

The damping of revivals can be partially explained by the expected variation in U2 over a non-uniform lattice due to an 
additional harmonic term in the trapping potential. Inhomogeneity in U2 causes dephasing due to the variation in the revival 
times for coherent states at different sites, however, the estimated 3-5% inhomogeneity of U2 should allow as many as 10-20 
revivals compared to the ~5 seen in ||5] |6l I?) . In contrast, we show below that coherent three-body interactions can cause 
dephasing of coherent states at each lattice site after only a few revivals. 

The effective theory in this paper describes the low-energy, small scattering length, small atom number per lattice site regime, 
for deep 3D lattices with negligible tunneling. These approximations are reasonable for the experiments in ||5] |6l |2l . Extensions 
of the analysis might include tunneling, including second-order f29] and interaction driven |30| tunneling, and the incorporation 
of intrinsic higher-body interactions. Effective field theory has also proven to be an important tool in the large scattering length 
limit |12|. It would be particularly interesting to simulate the controlled breakdown of the effective theory developed here by 
increasing the scattering length or atom number, or by tuning other lattice parameters. Looking beyond the realm of atomic 
physics, our analysis suggests interesting possibilities for using optical lattices to test important mechanisms in effective field 
theory fJT]. 

In Section 2, we construct a multimode Hamiltonian H which we use to obtain an effective single-mode Hamiltonian Hgff. 
In Section 3, we describe the physical processes that generate higher-body interactions. In Section 4, we estimate the effective 
three-body energy. In Section 5, we show how the coherent three-body interactions modify the collapse and revival dynamics. 
Finally, we summarize our results in Section 6. 

II. EFFECTIVE THREE-BODY MODEL FOR NEUTRAL BOSONS IN AN OPTICAL LATTICE 

A many-body Hamiltonian for mass nia neutral bosons in a single spin state can be written as 

I ij^HQiidv+ \ /v>(r)^A^(r')y2(r,r')i^(r)^A(r')drdr' (1) 
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+ IJ^^ (r) (r') V'^ (r") V3 (r, r', r") (r) (r') (r") drdr'dr" + 

where Vm are intrinsic m-body interaction potentials, and Ha is the Hamiltonian for a single particle in the optical lattice. We 
set Vm>2 = to focus on the physics of effective interactions induced by V2. In experiments, the effect of intrinsic and effective 
interactions are both present. 

It is our goal to construct a low energy, effective Hamiltonian iJgff for describing a small number of atoms in the vibrational 
ground state of a lattice site, while incorporating leading-order corrections from virtual excitation to higher bands. In the 
quantum mechanical approach, Huang et al. |32| have shown that a local regularized delta-function potential V2 
^(3) (j _ J.'") (djdg) g, where g = |r — r'| , can be used to obtain the low-energy scattering for two particles. To go beyond 
the two-particle case, we find it convenient to instead use the renormalization methods of quantum field theory and the non- 
regularized delta-function potential 

l/2(r,r') = <725('^(r-r'). (2) 

We regularize the theory in perturbation theory by using a high-energy cutoff A in the sum over intermediate states, which is 
equivalent to using a regularized (non-singular) potential. We view A as a physical threshold beyond which the low-energy 
theory fails. We note that the low-energy physics does not, in the end, depend on the method of regularization, and that the 
physical results found below after renormalization are insensitive to A. The key observation is that even if a fully regularized 
form of V2 is used renormalization is still required recognizing that the bare parameter 52 is not the physical (renormalized) 
coupling strength 52- (In the following we use a tilda to distinguish between bare and renormalized parameters.) 
Employing renormalized perturbation theory ||3T1 . we write 52 = 32 + c, where 



52-^^^^ + 0(aL), (3) 



is chosen to reproduce the exact, low-energy limit given in 11331 for two atoms in a spherically symmetric harmonic trap, and 
flscat is the scattering length at zero-collisional energy. The first-order approximation to 32 suffices for the calculation of the 
three-body energy at second order given below. The value of the counter-term c, which cancels the contributions to the two- 
body interaction energy that diverge with A, is determined by the normahzation condition Eq. (|3]l. The local Hamiltonian with 
counter-term and physical coupling parameter becomes 

f i^^Ho^jdv + ^ (52 + c) / i^^^^dv. (4) 



3 



To develop a low-energy effective theory for a deep optical lattice, we expand the field over a set of bosonic annihilation 
operators a^^ and single particle wavefunctions (r) giving (r) — J2if_i '^v '^V' where the indices — {fi^, f^y, f^'z} 
with iix.y.z — 0, 1, 2, ... label 3D vibrational states and i labels the lattice sites. To focus on the role of interactions we assume 
a deep lattice with Ug > "i states per spatial dimension at each site, making tunneling of atoms in the ground vibrational state 
/i = {0, 0, 0} = negligible on the timescale of interest [5J. Since we are not considering the role of tunneling, for simplicity 
we use isotropic harmonic oscillator wavefunctions at each site with frequency lo and length scale a = ^JWfrn^ determined by 
the (approximately) harmonic confinement within a single lattice well. Note that even with tunneling neglected, anharmonicity 
of the lattice potential is a potentially significant effect. We also expect our model to break down or to require significant 
modification for very shallow lattices or near the Mott-insulator phase transition where the effects of tunneling are important 
Il2i[34l[35l|36l. 

Inserting the expansion for into Ji, interchanging the order of integration over r and summation over modes, and dropping 
terms that transfer atoms between sites (e.g. tunneling), we obtain for each lattice site the multimode Hamiltonian H = Hq + H2 , 
where 

Ho = Yl (5) 



H2^ \ (U2 + A)y^ K^i,c,\a^nO'la„ax. (6) 



and 

Ho ^- (u-, + A\Y 

For brevity we suppress the lattice site index i. The single particle energies are = {fi^ + l-hj + f^z) setting the ground 
state energy = i?{o,o,o} = 0- The two-body interaction energy for ground state atoms is 

U2 = %3 — = \ -huj (ascat/o-) , (7) 

and A = (27r)^'^^^ c/a'^ is the counter-term in units of energy. The matrix elements 

K^uiS = (27r)^^^ (T^ / (t)f,(t)^4>^(j)sdr (8) 



are normalized so that Xoooo = and they vanish for transitions that do not conserve parity. It should be noted that, when 
there is a cutoff in the sum over modes, both the regularized and non-regularized delta-function potential lead to the same 
Hamiltonian H and matrix elements in Eq.ljsjl, and thus they produce the same results in the regularized (cutoff) quantum field 
theory. We emphasize that after the renormalization of the two-body interaction energy, the induced three-body interaction 
energy is insensitive to the cutoff A. We develop the perturbation theory in the small parameter ^ defined by 

U2 / 2 flscat 

huj V TT (7 



e^^-Vz— +OKat)- (9) 



The total interaction energy for n atoms in the vibrational ground state in the single mode per site approximation is Ei^t = 
U2n {n~ 1) /2. Commonly, a single-mode approximation is made based on the two-body interaction energy per particle being 
much less than the band gap, i.e., E^^^/n = {/2 (n — 1) /2 <C fiLjj or ^ 1. For ^^Rb with scattering length Oscat — 5.3 run 
and a lattice with ijj/2ti ~ 30 kHz, we have [/2//1 ~ 2.0 kHz, and ^ = 0.07. We will use these as typical system parameters in 
the following analysis. With ^ — 0.07, the single mode per site condition n ^ ~ 15 is easily satisfied and the influence of 
higher-bands will produce only small (though important) corrections. For coherent states, for example, we show that small three- 
body energies can lead to large phase shifts over time resulting in interferometric-like sensitivity to higher-body and higher-band 
processes. 

To obtain an effective Hamiltonian iJeft, we use the multi-mode Hamiltonian H = Hq + H2 to compute the atom-number 
dependent energy shift for atoms in the vibrational ground state. Our approach is essentially equivalent to the effective field 
theory procedure of summing up to a cutoff over all 'high-energy' modes p. with > /lo;, which generates a low energy 
effective theory with all consistent local interactions. We obtain an effective Hamiltonian ifeft for the /i = mode that is vaUd 
in the low-energy regime Ei^Jn ^ nU2 ^ fuxi, which is consistent with the single mode approximation discussed above. Of 
course the multimode Hamiltonian H itself is an effective Hamiltonian which is only valid for energy scales E^ + Eint/n <C 

In the case of atoms confined in a deep well, the effective Hamiltonian for U2 ^ huj is 

H,s = Eoa^a + ^ i/^a^^a"/™!, (10) 

m>l 



4 




FIG. 1: The effective two-body interaction energy U2 is given through second order by diagrams (a)-(d). Diagram (d) is the counter-term 
that cancels the diagrams (b) and (c), fixing 172 as the physical (renormalized) two-body energy. Diagrams (f)-(i) are examples of processes 
contributing to the effective three-body interaction energy f/3, represented by diagram (j). Diagram (g) gives the leading order contribution, 
assuming Us = 0; it shows how an effective three-body interaction involving three distinct incoming particles arises at second order in 
perturbation theory. Diagrams (h) and (i) are two of the effective three-body processes that arise at third order (others are not shown). If the 
bare three-body vertex shown in (f) does not vanish additional three (and higher) body counter-terms are also required. 



where creates an atom in a renormalized ground vibrational state. The Eqcl'o. term vanishes since we set £'0 = 0. The 
dominant term in i/eff is the two-body energy, and the higher-body interaction energies scale as nU„i/tj„i_i ^ {nU2/fu^i^ ^ 
< < 1. ^ 

The energies [/,„ can be computed in perturbation theory in the small parameter ^ using H to find the energy of n atoms 
in the ground vibrational mode. At m}^ order in ^, all local interactions up through the {m + l)-body term i/,n+i = 
Um+i(i^™'^^oJ"^^ / {in + 1)1 are generated. In this paper, we work to second order in ^ for which the effective Hamiltonian 
is 

Heff = U2a^^a^/2 + U^iO^^a^ /Q. (11) 

Using fi = a^a and [a,a^] = 1, the two- and three-body terms can be written as a^'^a'^ = n{n—l) and a^^a^ — 
n(n — 1) (n — 2) ; the latter expression shows explicitly that the effective three-body interaction only arises when there are 
three or more atoms in a well. Eigenstates of ifeff with n atoms have energies 

E (n) = U2n {n - 1) /2 + U^n {n - 1) {n - 2) /6. (12) 

Note that the three-body energy scales as and thus its influence relative to the two-body term, though small, can be tuned by 
changing the number of atoms in a well. 



III. MECHANISM FOR EFFECTIVE INTERACTIONS 



We now describe the virtual processes that give rise to effective m-body interactions in a deep lattice. Writing the perturbative 
expansion for the energy of an n atom state |n) through second order as £^ (n) — i?'"^ {n)+E^^'> + (n) , the zeroth-order 
energy is i?^") (n) — Equ — 0, recalling that Eq — 0. The first-order energy shift, treating H2 as the perturbation Hamiltonian, 
is the usual expression 

£;(!) (n) = (n\ H2 \n) = U2n (n - 1) /2. (13) 

This is the leading order result for the two-body interaction energy and, setting n — 2, the renormalization condition Eq. (j?]) 
shows that yl = to first order in U2- Figure 1(a) represents this first-order process. 
The second-order energy shift can be written as 

(n) = X: +Anin- 1) /2 (14) 

with K^i, = K^yQQ and 11 > v. The O (ui^ counter-term A now appears. At this order, A is determined by the renormalization 
condition £'(2) = £(°) (2) + (2) + £(2) (2) = U^, implying that £(2) (2) = 0. The sum is over intermediate states 
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l/ii^) = ajj^alaouo \n) with energy E/^i, = + > this excludes the = v = Q state. For regularization purposes we 
introduce a high-energy cutoff that limits the sum to E^y < A. The factor s^,y = {4, 1} if {fi = ly, fj, ^ comes from the 
two equivalent terms aj^aj, and ala)^ that appear in H2 - Each term in E^'^'^ involves a two-body collision-induced transition to a 
virtual intermediate state. For example, the state llo^lx) corresponds to two atoms both excited along the x direction with energy 
Ell — (note that Kii — 1/4 for this transition), with the remaining n — 2 atoms left in the /i = mode. Because collisions 
conserve parity, contributions from states like jl^^lj,) vanish. 



The crucial observation is that the series in Eq. ( 14 1 separates into two distinct sums corresponding to two-body and three-body 
interactions, respectively, i.e., 

£;(2) {n) = 6U2n {n - 1) /2 + SU^n {n - 1) (n - 2) /6, (15) 



where SU2 includes the counter-term contribution A from Eq. ( [T4| . For fi ^ and v ^ Q intermediate states, 
I (/ii^l aJ^aJ,aoao |n) | = kti (71 — 1) where n — {2, 1} if {/i — t^, 7^ i/} , with the factor of 2 resulting from Bose stimula- 
tion when both atoms transition to the same excited state. Because these terms are proportional to n (n — 1) they contribute to 
the two-body energy shift 5U2- A diagram representing this two-body process, with two atoms colliding, making transitions 
to virtual excited vibrational states, and then returning to the ground state after a second collision with each other, is shown in 
Fig. 1(b). The fJ, virtual states and /i = vibrational ground states are represented by dashed and solid lines, respectively. 
The origin of the three-body energy can be seen by examining the ^ > 0,1/ — intermediate states. We have 

I (/iz^l aj^ajaoflo \n) — n{n — 1)^ = n {n — 1) + n {n — 1) {n — 2) , (16) 



showing that these terms generate both effective two- and three-body energies. The extra factor of (n — 1) in Eq. ( 16 1 results 



from Bose stimulation of an atom back into the /i = state when two atoms collide but only one makes a transition to an excited 



state. Figure 1(c) shows the two-body process corresponding to the n (n — 1) term in Eq. ( 16 1. Figure 1(d) shows the counter- 
term A whose value is determined such that it cancels the contributions from Figs. 1(b) and (c), thereby maintaining, through 
second order, the renormalization condition that the parameter U2 is equal to the physical two-body energy. To arbitrary order 
the renormalization condition determines A such that all higher-order two-body diagrams cancel, as represented by Fig. 1(e). 



Figure 1(g) shows the effective three-body process corresponding to the n (n — 1) (n — 2) term in Eq. ( 16 1. This process 
gives the leading-order contribution to SU3 and generates a three-body interaction energy C/3 = t/3 + 6U3 even if the bare U3, 
represented by Fig. 1(f), vanishes. More generally, we expect U3 ^ 0, but nevertheless the contribution to given by 8Uz can 
be a significant (possibly even dominant) correction. Looking at Fig. 1(g), we see that two initial /i = atoms collide giving 
rise to one n ^ atom that subsequently collides with a third, distinct fi = atom. In Fig. 1(g) there are three distinct incoming 
atoms resulting in an effective three-body interaction mediated by the fJ, ^ intermediate state. The renormalized three-body 
interaction energy is represented in Fig. l(j) by a square vertex with three incoming and outgoing particles. Figures 1(h) and 
(i) show examples of two different processes contributing to U3 at third-order in ^; they illustrate how higher-order processes, 
including counter-terms, arise. Their contributions, and other third-order processes not shown, are not explicitly computed 
below. At third order, effective four-body interactions also arise. 

Notice that there are two types of diagrams in Fig. 1: tree diagrams [e.g. Fig. 1(g)] and loop diagrams [e.g. Fig. 1(b)]. In 
general in quantum field theory the contributions from some loop diagrams diverge with the cutoff A, necessitating the need for 
renormalization, whereas the contributions from tree diagrams are finite 1 3 1 1 . We will see this behavior explicitly below. In fact, 
at m*'' order in ^, there will be a set of tree diagrams giving a. finite, leading-order contribution to the effective (m + l)-body 
interaction energies Um+i- We note that even if all intrinsic higher-body interactions exactly vanish there will be effective m- 
body interactions and associated energy scales [/„i generated by the two-body interactions. Consequently, the nonequilibrium 
dynamics of n atoms in the ground vibrational mode, when ^ 1, will be characterized by a hierarchy of frequencies 

(U2/h,U3/h...,Um/h] 



IV. ESTIMATE OF THE EFFECTIVE THREE-BODY INTERACTION ENERGY 



Returning to Eq. ( 14 1 for the second-order energy shift and separating it into two- and three-body parts, we find that 

SU2^~U^(Y!' KlJE^^+A, (17) 



and 



6U, = -6Ui (Y^l^^ KlJE^o) ■ (18) 
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In the expression for SU2 the sum is over all fi and ly (both ii> v and v > ji) except for the ^ = v = Q mode. Similarly, in the 
expression for SU^ all ^ except for /i = are summed over. 

As expected, the sum X^^iy ^fny/^p-f corresponding to the second order, 1-loop diagram in Fig. 1(b) diverges with A, re- 
flecting the divergent relationship between the bare U2 and renormalized U2 energy parameters. In fact, the sum scales with the 
cutoff as A^/^. The renormalization condition that E (2) = U2 determines A by requiring that SU2 = 0. To second-order, the 
interaction energy of n atoms is thus 

E (n) = U2n {n - 1) /2 + SUsn {n ~ I) {n - 2) /6, (19) 

assuming U3 = SU3. 

After cancelling the two-body corrections with A, the remaining second-order term gives an induced three-body energy that is 
insensitive to A : the quantity Xl^>o ^^o/^mO corresponding to the second-order tree diagram in Fig. 1(g) converges. Writing 

SUs/hio = -/3e^ (20) 
this sum can be solved analytically for a spherically symmetric harmonic trap in the A ^ 00 limit, and we find IIJTII 

/3 = 4\/3-6 + 61og f — ^ I -1.34.... (21) 

Cutting off the sum at E^^/hio < A/huj = 4 akeady gives f3 ~ 1.30 showing the rapid convergence of the series. The 
convergence of this sum is an example of the generic behavior that contributions from tree diagrams are finite. If the bare is 
zero or sufficiently small, the effective three-body energy is negative, giving attractive three-body interactions, and reducing the 
total interaction energy for both positive or negative U2- 

We expect significant corrections due to the anharmonicity of the true lattice potential. The single-particle energies of higher 
vibrational states are lowered on the order of the recoil energy Eji, defined as the gain in kinetic energy for an atom at rest that 



emits a lattice photon. This leads to a decrease of the energy denominator in Eq. ( 14 1 and, for the typical system parameters 
considered here, this can give an estimated correction to C/3 of 10% or more. The matrix elements K^i, will also have corrections. 
These effects can be computed numerically using single-particle band theory. 

We have defined our perturbation theory around the zero-collisional energy limit, but in a trap the collision energy of ground 
state atoms is on the order of huj. As shown in |38| , an improved treatment replaces the zero-energy scattering length by an 
effective scattering length defined as 

- — = — + lr,k^, (22) 

where the effective range Tg is on the order of the van der Waals length scale (maCe/^i^)^^^ away from a Feshbach resonance, 
and the collision energy is fi^fc^/ma For *^Rb the van der Waals length is approximately 8 nm. In a trap the ground vibra- 
tional state wavevector k ~ cr^^ produces a fractional increase in scattering length on the order of {re/cr) C By incorporating 
the effective scattering length model we can extend the range of validity of our model. 

Even neglecting these corrections, the perturbation theory generated by Eqs. Q and ^ does not predict the two-body energy 
U2 but instead uses the measured value, or the exact result calculated by other methods such as Busch et al ||33l . as input from 
which SU3 is obtained. Similarly, the effective theory does not yield the intrinsic three-body interaction energy U3, and therefore 
^^3 = ^3 + must also be determined by either measurement or a theory of three-body physics if the intrinsic interaction 
energies Um>2 are non-zero. On the other hand, the effective theory shows that even if ?7m>2 = there are significant induced 
three- and higher-body interactions, and if non-zero f7„i>2 measured the effective contribution from two-body processes 
must be taken into account before the intrinsic higher-body coupling strengths can be extracted. Note that if non-zero bare 
(intrinsic) parameters Um>2 are included in our model, additional counter- terms will be needed to cancel divergences, reflecting 
the need to ultimately determine any intrinsic higher-body coupling strengths via either measurement or an exact high-energy 
theory. 

Assuming Us ~ 0, Fig. 2 shows U2 = (.fiuj and versus ^, including positive (^ > 0) and negative (f < 0) scattering 
lengths. Using ^ ~ 0.07 for ^^Rb in a 30 kHz well gives U2/h ~ 1.9 kHz and Uj,/h ~ —200 Hz. Using a Feshbach resonance 
ll40l to change Oscat and thus ^, or fixing Oscat and changing the trap frequency w, it is possible to tune the relative strengths of 
the three-body (and higher-body) interactions. It would be interesting to explore the breakdown of the perturbative model by 
increasing either ^ or the atom number n, or by decreasing the lattice depth so that the influence of tunneling and higher-band 
effects increases. 



7 




FIG. 2: The figure shows Us and U2, in units of hcj, versus ^. The bold hne shows the induced three-body energy SUs = U3 = -13^'^ with 
/3 = 1.34, assuming the intrinsic energy Us vanishes. The dashed line shows the leading-order two-body energy U^- The graphs extends 
beyond the regime of strict validity of the perturbation theory, which requires ^ 1 where n is the number of atoms in a lattice well, to 
illustrate the overall scaling of the two- and three-body energies. The collapse and revival experiments in |5]|6||3 have uj/2n ~ 30 kHz and 
£, ~ 0.07, putting them well within the perturbative regime. The inset shows U3 for the range < ^ < 0.1. 




FIG. 3: Collapse and revival visibility versus time t, with ^ = 0.07 and n = 2.5. Curve (i) shows the case with neither inhomogeneities nor 
three-body interactions included. Curve (ii) shows the effects of ~ 5% inhomogeneities in U2- Curve (iii) shows the effects of three-body 
interactions with P — 1.34 but no inhomogeneities. Note that the three-body mechanism influences the visibility of revivals immediately, and 
it will be important even if inhomogeneities are stronger than are shown in curve (ii). Curve (iv) shows the combined effects of inhomogeneities 
and three-body interactions. 



V. DYNAMICS AND DECOHERENCE OF ATOM-NUMBER COHERENT STATES 



We now investigate the influence of effective three-body interactions on the phase coherence of an N atom nonequilibrium 
state (0)) — (X]f=i ^lo 1^^) /VM)®'^, obtained by quickly loading a BEC into a lattice with M sites. To a good approxima- 
tion the state can be treated as the product of coherent states lISl fTOl . 

i i 

where \ai) — Ui \ai) and \ai\^ — fii is the average number of atoms in the i^^ site. A relative phase (pij between sites 
i ^ j exists when {a\aj) = rje^'^'^ and rj ^ 0. The initial state |^ (0)) has ij — n, and there are well-defined relative phases 
= for all i,j in this case). In contrast, the equilibrium Mott insulator state, achieved by much slower loading f35', 'SSI, 
has approximate number states in each well giving 77 « 0, though there can be some degree of short-range phase coherence 

imiEiisi. 

Coherent states in optical lattices make natural probes of higher-body coherent dynamics because small atom-number de- 
pendent energies can lead to significant phase shifts over time. After a hold time th in the lattice, the initial state evolves to 
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I* {th)) ~ III 1^7 ith))i , where the state of the i*^ well is 



„ V n! 



and Ei (n, ) is given in Eq. ( 12 1, restoring the index i labeling the lattice site. Snapping the lattice off at time th, the wavefunctions 
from each well freely expand for a time until they fully overlap, analogous to the diffraction of light through a many-slit 
grating. 

For a uniform lattice, the fringe visibility is ||5l 



V{th)^\{v{th)\a\v{th))f /n. (25) 

With no inhomogeneities and setting = 0, we obtain V (th) = e-'^"-[^-''°''{^2*h/ti)] j^ie visibihty for n = 2.5 is plotted as 
the thin dashed line labeled (i) in Fig. 3, showing the well-known collapse and revival dynamics with period <2 = h/U2. For the 
^''Rb system parameters used here t2 — 0.52 ms. 

The thin line labeled (ii) in Fig. (3) shows the influence of an approximate 5% variation in the two-body energy U2- We 
average the (t) over a 60 lattice-site diameter spherical distribution. While the effect of inhomogeneities are important, a 
larger variation in U2 then expected would be required to explain the decay of interference fringes after only 5 revivals as seen 
in experiments ||5]|6]|7]. We note that the longer timescale for three -body recombination can be distinguished from the coherent, 
number conserving interactions derived here by tracking changes in total atom number, and this appears to be negligible on the 
revival damping timescale |26 1. Other mechanisms, such as non-adiabatic loading |44 45 1 and collisions during expansion B6l 
will reduce the initial fringe visibility but do not explain the rapid decay of the visibility versus hold time th. 

To compute the visibility with three-body interactions we numerically evaluate 

(t,)\a Iry (t,)) = ae'^ V ^,-H^.+u,i-^i)/2]t./n ^ ^26) 



n=0 



The bold (blue) dashed line labeled (iii) in Fig. 3 shows the visibility V (th) = |(?7 {th)\ a, \ri {th))\^ /n versus tfi/t2 assuming 
no inhomogeneities, n = 2.5, and the harmonic oscillator value (3 = 1.34.... With ^ = 0.07, C/3 = 0, and llj/2tt = 30 kHz, the 
effective three-body frequency is Us/h ~ —200 Hz, and 1/2/^ — 2.1 kHz. The relatively small effective three-body interactions 
have a strong effect on the coherence of the state and the resulting quantum interference, showing how collapse and revival 
measurements can be a sensitive probe of coherent higher-body effects. The dephasing is faster than may have been expected 
from the small size of U3 because the three-body energies scale as Usn^ versus [72^^ for two-body energies, and thus have an 
increased influence on higher-number components of a coherent state. Similarly, coherent states with significant n > 4 atom 
number components will probe the four- and higher-body interaction energies. The bold (red) solid line labeled (iv) in Fig. 3 
shows the combined effect of both ^5% inhomogeneities in U2 and three-body interactions. 

The decay of the visibility in Fig. (3) is faster than what is seen in |5 , 6, 7 1. Figure 4 illustrates the sensitivity of the evolution 
of the visibility to the three-body energy scale by showing three cases corresponding to U^/h = {—200, —150, —100} Hz. The 
curves have been displaced vertically for clarity. Curve (i) for Us/h — —200 Hz corresponds to U3 — 0, f3 = 1.34, ^ = 0.07, 
and w/27r = 30 kHz. Curve (ii) corresponds to a reduced U^/h = —150 Hz, which could be the result, for example, of a positive 
intrinsic three-body energy Us/h — 50 kHz, or a change in parameters giving either /3 ^ 3/3/4 or ^ ^ \/3C/2- Similarly, curve 
(iii) corresponds to = 100 Hz, which could be due to a positive intrinsic three-body energy U^/h — 100 kHz, or to a change 
in parameters giving either /? — > /3/2 or ^ — > ^/\/2. The collapse and revival visibilities are also very sensitive to the average 
atom number n. A smaller value of {73 appears to agree better with the initial damping seen in 15] IS 13, and this may indicate 
the presence of a non-zero intrinsic U3. However, accurate measurement of the system parameters is necessary if a value of 
the intrinsic U3 is to be obtained using C/3 = [/g — 6U3. Nevertheless, it is clear from Fig. (4) that both intrinsic and induced 
three-body interactions can be important on experimentally relevant timescales. 

Figure 4 also shows the partial and full revivals resulting from the beating between two- and three-body frequency scales 
expected if inhomogeneities are sufficiently reduced. The period for nearly full three-body revivals <3 = h/Us gives a direct 
method of measuring U3. Recently, long sequences of collapse and revivals showing multiple frequencies have been reported 
li26i ; our analysis suggests that these may be used to study higher-body interactions in optical lattices. 



VI. SUMMARY 



We have shown that two-body induced virtual excitations of bosons to higher bands in a deep 3D optical lattice generate 
effective three-body and higher-body interactions. Although our methods do not yield the intrinsic higher-body interaction 
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Two-body(i) = 200 Hz (e.g., = 0) 
revivals 



'c/2 



Three-body 
revivals 




(ii) [7,= 150 Hz (e.g., (7, = 50 Hz) 





(iii) f75=100Hz (e.g., (7, = 100 Hz) 




FIG. 4: The figure shows the collapse and revival visibility versus time for Us = 200, 150, and 50 Hz assuming negligible inhomogeneities. 
Curve (i) corresponds to [/s = 0,^ = 0.07, /3 = 1.34, n = 2.5, and ui/2-k — 30 kHz. Curves (ii) and (iii) correspond to smaller three- 
body energies U3, which could be due to a non-zero intrinsic U3, a reduction in /3, or a change in other system parameters including ^ or n. 
Three-body revivals occur at multiples of = h/Us, providing a method for measuring the coherent three-body interaction energy. 



energies Um, we find that even if U„i — there are significant effective interactions that can have a surprisingly strong influence 
on the dynamics of non-equiUbrium coherent states. The mechanism for higher-body interactions is based upon the recognition 
that at low energies the presence of excited (i.e. higher-energy) vibrational states manifest as m-body terms in an effective 
Hamiltonian iJeff • While it is possible for an effective (or renormalized) m-body interaction to vanish or to be very small due 
to close cancellation of the intrinsic (i.e. Um) and induced (i.e. SUm) contributions to Um, we do not expect this to happen 
in general. It is possible to tune the relative effective m-body interactions by exploiting Feshbach resonances to control ascat> 
or by changing the lattice potential. This suggests intriguing possibilities for probing and controlling the physics of effective 
field theories (e.g., effective interactions, running coupling constants, and the emergence of non-perturbative effects) in optical 
lattices. Using optical lattices to simulate the controlled breakdown of an effective field theory would be particularly interesting. 

We would like to thank J. Sebby-Strabley and WD. Phillips for very helpful conversations. P.R.J, acknowledges support from 
the Research Corporation for Science Advancement. J. V.P. acknowledges support from lARPA. 
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